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Exploring accretion theory with X-ray binaries in the SMC 
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ABSTRACT 

The understanding of the accretion process on to compact objects in binary sys- 
tems is an important part of modern astrophysics. Theoretical work, primarily that 
of Ghosh & Lamb (1979), has made clear predictions for the behaviour of such sys- 
tems which have been generally supported by observational results of considerably 
varying quality from galactic accreting pulsar systems. In this work a much larger 
homogeneous population of such objects in the Small Magellanic Cloud (SMC) is used 
to provide more demanding tests of the accretion theory. The results are extremely 
supportive of the theoretical models and provide useful statistical insights into the 
manner in which accreting pulsars behave and evolve. 
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1 INTRODUCTION AND BACKGROUND 

The Be/X-ray systems represent the largest sub-class of all 
High Mass X-ray Binaries (HMXB). A survey of the litera- 
ture reveals that of the ~240 HMXBs known in our Galaxy 
and the Magellanic Clouds (Liu et al., 2005, 2006), >50% 
fall within this class of binary. In fact, in recent years it has 
emerged that there is a substantial population of HMXBs in 
the SMC comparable in number to the Galactic population. 
Though unlike the Galactic population, all except one of the 
SMC HMXBs are Be star systems. In these systems the orbit 
of the Be star and the compact object, presumably a neu- 
tron star, is generally wide and eccentric. X-ray outbursts 
are normally associated with the passage of the neutron star 
close to the circumstellar disk (Okazaki & Negueruela 2001), 
and generally are classified as Types I or II (Stella, White 
& Rosner, 1986). The Type I outbursts occur periodically 
at the time of the periastron passage of the neutron star, 
whereas Type II outbursts are much more extensive and oc- 
cur when the circumstellar material expands to fill most, 
or all of the orbit. This paper concerns itself with Type I 
outbursts. General reviews of such HMXB systems may be 
found in Negueruela (1998), Corbet et al. (2008) and Coe et 
al. (2000, 2008). 

This paper reports data acquired over 10 years using the 
Rossi X-ray Timing Explorer (RXTE) on the HMXB popu- 
lation of the SMC. During the period of these observations 
there have been many opportunities to study spin changes 
arising from accretion torques. This extremely homogeneous 
population permits the first high quality tests to be carried 
out of the work of Ghosh & Lamb (1979) and Joss & Rappa- 
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Figure 1. An example of the spin period changes seen in one 
system, SXP59.0. The straight lines show the best fit to each 
data set. 



port (1984). The simplified source naming convention used 
in this work follows that established by Coe et al (2005). 



2 OBSERVATIONS 

The SMC has been the subject of extensive monitoring using 
the RXTE Proportional Counter Array (PCA) over the last 
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10 years. The PCA instrument has a Full Width Half Max- 
imum (FWHM) field of view of 1 ° and data in the energy 
range 3-10 keV were used. Most of the observations were 
pointed at the Bar region of the SMC where the major- 
ity of the known X-ray pulsar systems axe located. Sources 
were identified from their pulse periods using Lomb-Scargle 
(Lomb 1976, Scargle 1982) power spectral analysis of the 
data sets. Laycock et al. (2005) and Galache et al. (2008) 
present full details of the data analysis approach that has 
been used to determine which pulsars were active during 
each observation. In their work, for each X-ray outburst, 
the pulse amplitude and history of the pulse periods were 
determined. Those results are used in this work. Since a 
database of >10 years of observations of the SMC exists it 
was therefore possible to use these data to search for evi- 
dence of spin period changes in the systems. The PCA is 
a coUimated instrument, therefore interpreting the strength 
and significance of the signal depends upon the position of 
the source within the field of view. In all the objects pre- 
sented here the target was assumed to be located at the 
position of the known optical counterpart. Only observa- 
tions that had a collimator response >25% and a detection 
significance of >99% were used in this work. 

A total of 24 sources were chosen for this study. In each 
case three possible measurements of period changes were 
obtained: 

• Individual active periods lasting typically 50-500 days 
were studied and used to determine the P for a particular 
source data subset (referred to in this work as the ShortP). A 
simple best fit straight line to the pulse history plot was de- 
termined using a minimising technique. No attempt was 
made to fit more complicated profiles to the data, though in 
some cases higher order changes are suggested. An excellent 
example of the spin period changes seen in these systems 
is presented in Figure [T] which shows two outbursts from 
SXP59.0. Clearly both outbursts indicate an initially higher 
P which levels off as the activity period progresses, but only 
the weighted average is used in this work. One other fac- 
tor that could also modify the spin period history would 
be Doppler-related changes. However, attempts to fit period 
histories with binary models have always proved difficult 
(but see, for example, Schurch et al, 2008, for one possi- 
ble success) suggesting that the changes are dominated by 
accretion driven variability. 

• In addition to the ShortP, where possible a longer term 
value was determined for each source from the whole data 
set covering ~10 years of observing - referred to here as 
the LongP. This typically included several periods of source 
activity with significant inactive gaps in between. 

• For many sources the orbital period is clearly apparent 
in the sequence of X-ray outbursts. For others, optical data 
from the Optical Gravitational Lensing Experiment (OGLE) 
project (Udalski, Kubiak and Szymanski 1997) have been 
used with good success to discover the orbital modulation - 
see Section 3.4 below for further discussion on this point. 

It was not always possible to determine both a ShortP 
and LongP for every source in this work due to several 
possible reasons; one being the observational coverage and 
another the activity history of the system. Details of the 
recorded spin period changes are given in Table [T] Full 
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Figure 2. The observed relationship between equilibrium spin 
period and the rate of spin change during an outburst (Short P). 
The straight line shows a P = kP^ relationship. The four points 
marked with a solid squares have positive P values (spin down), 
the rest are negative (spin up). 



records of the behaviour of each source may be found in 
Galache et al. (2008). 

The strong link between the equilibrium spin period and 
the rate of change of spin period seen during outbursts is 
shown in Figure[2] In this figure, the straight line represents 
P — kP^ - as predicted for accretion from a disk on to a 
neutron star (see Equation 15 in Ghosh & Lamb, 1979). 
It is interesting to note that the four spin-down systems 
(SXP8.80, SXP59.0, SXP144 & SXP1323) fit comfortably 
on this relationship alongside the much larger number of 
spin-up systems. 

To pursue the understanding of the accretion process 
further, a value for the X-ray luminosity, , is needed dur- 
ing each outburst. So the X-ray luminosity was calculated 
from the observed peak pulse amplitude in counts/pcu/s 
that occurred during the outburst that produced the ShortP 
values listed in Table [T] This amplitude is converted to 
luminosity assuming that 1 RXTE count/pcu/s = 0.4x10'^^ 
erg/s at the distance of 60kpc to the SMC (though the 
depth of the sources within the SMC is unknown and may 
affect this distance by up to zblOkpc). The X-ray spectrum 
was assumed to be a power law with a photon index — 1.5 
and an Nh = 6 x lO'^" cm~^. Furthermore it was assumed 
that there was an average pulse fraction of 33% for all the 
systems and hence the correct total ffux is 3 times the 
pulse component. Thus the values shown in Table [T] were 
determined using the relationship: 

Lx = 0.4 X 10^'' X 3P erg s"^ 

where R = RXTE counts PCU"^ s"^ 

Note that though the values for obviously scale lin- 
early with pulsed fraction, the neutron star mass determi- 
nations discussed below in Section 3.3 are very insensitive 
to any value between 10% and 60% for the pulsed fraction. 

Figure |3] (upper panel) shows a histogram of all the 
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Table 1. Pulse period determinations. The X-ray luminosity (3-10 keV)comes from the peak amplitude of the pulsed modulation signal 
during an outburst - see text for details. The pulse period given in column 2 is the average value during the periods of activity used for 
measuring Short P. Two systems (SXP59.0 and SXP172) have had two separate outbursts measured and so appear twice in the table. 
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known pulse periods for accreting pulsar systems in the SMC 
(see Galache et al, 2008 for details of the vast majority of the 
sources). Because of the length of the each observation (typ- 
ically lOks - 15ks) there is undoubtedly an instrumental cut- 
off starting around >lks affecting the ability to detect longer 
periods. The sharpness of this cutoff is partially a function 
of the pulse profile. At the short period end, the data were 
regularly searched down to periods of 0.5s (Galache et al, 
2008) and should be complete down to that number. The 
lower panel in Figure |3] shows the distribution of known bi- 
nary periods determined from either the sequence of X-ray 
outbursts and/or a coherent modulation of the optical hght 
from the system. 

Figure 0] shows a histogram of the peak outbursts ob- 
served from the systems studied in this work. The signifi- 
cantly larger outburst seen from SXP2.37 is excluded from 
this plot - this one outburst may well be tending towards a 
Type II outburst. 



3 DISCUSSION 
3.1 Reverse torques 

It is interesting to note that of the 15 systems for which 
short term P changes were measured, 4 (SXP8.80, SXP144, 
SXP59.0 and SXP1323) show spin-downs during outburst 
rather than spin-ups. Such reverse torques have been re- 
ported before and even transitions from spin-up to spin- 
down observed for example, 4U1626-67 (Bildsten et al, 
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Figure 3. The histograms of the observed pulse periods (upper) 
and binary periods (lower) for SMC systems. The majority of 
these systems are described in Galache et al. (2008), but others 
from the published literature have also been included. 
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Figure 4. The histogram of the estimated peak X-ray luminosi- 
ties for the outbursts discussed in this work. The exceptionally 
large outburst from SXP2.37 is excluded from this plot. 



1997). No such transitions are observed in the sample pre- 
sented here, but the numbers do give a clear indication of 
the frequency with which these reverse torques occur. In ad- 
dition, it is obviously a phenomenon that is not restricted to 
a particular spin period regime since our three object have 
periods ranging from 8.9s to 1325s. Five further systems 
show long term spin-down changes - see Table [1] However, 
all these other systems show spin-up during outbursts and 
the long-term change may be attributed to the gradually 
slowing down of the neutron star in the absence of repeated 
accretion episodes, rather than reverse torques. 



3.2 Spin period evolution 

The relationship between short spin period changes (seen 
during an outburst of typically 1-2 months) and the long 
term period changes seen over ~10 years worth of study is 
shown in Figure [S] The sources broadly follow a relationship 
indicating that LongP = 0.2 x ShortP. So individual out- 
bursts typically spin up (or down) a neutron star ~5 times 
faster than the longer, time-averaged changes. In reality, the 
neutron star is subjected to a series of "kicks" during each 
short outburst, followed by an intervening recovery period. 
This effect is illustrated clearly in Figure [T] 

3.3 Neutron star masses 

Ghosh & Lamb (1979) established the relationships for ac- 
cretion onto neutron stars formulating the now well known 
relationship between P and X-ray luminosity (see Equation 
15 in their paper). This relationship was re-iterated by Joss 
& Rappaport (1984) and others, with all sets of authors 
comparing the theoretical predictions with available data 
from galactic accreting pulsar systems. The quality of the 
observational data used was very variable, with many ob- 
jects represented simply by upper limits. In this work the 
consistently higher quality of the measurements allows us to 
compare theory and observation much more precisely. This 
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Figure 5. The relationship between the absolute values of the 
short (^couple of months) spin changes and the longer term 
changes (^10 years). All 11 sources from Table 1 with both short 
& long term P values are included. The dashed line shows the 
simple relationship LongP = 0.2 X ShortP. 



comparison is shown in Figure |6] where it is immediately 
obvious that the data provide strong support for the mod- 
els over two orders of magnitude. The probability of get- 
ting this level of correlation with an uncorrelated data set 
is 2.5x10"^ 

From the results presented here there is a strong sup- 
port for the average mass of neutron stars in these systems to 
lie between 1.3Mq and 1.9Mq. In fact, if the sample shown 
in Figure |B] is used, then the average mass of the neutron 
stars in the SMC sample is found to be 1.62±O.29M0. The 
only significant deviation from the average value is that of 
SXP7.78 (=SMC X-3) which strongly suggests a much lower 
neutron star mass for that system. SXP59.0 has the largest 
error bars of all the points and lies ~2(j away from the 
M=1.9Mq; hence is probably not significant measurement 
of a large neutron star mass. 

3.4 Implications for HMXB evolution 

There are currently 56 systems with known pulse periods, 
but only 27 with confirmed binary periods. Many of the bi- 
nary periods have been determined or confirmed from OGLE 
III long term data which frequently show evidence of mod- 
ulation believed to be at the binary period. Results such as 
those for SXP46.6 (McGowan et al, 2008) demonstrate the 
strength of this link by revealing regular optical outbursts 
coincident with X-ray outbursts. However, for many of the 
SMC systems the binary period has yet to be confirmed and, 
in particular, the longer term periods (>l-2 years) are in- 
creasingly difficult to ascertain. This may be due to either 
missing X-ray outbursts or/and a lack of optical modulation 
due, perhaps, to the size and shape of the orbit. 

However the period histogram should reflect accurately 
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the distribution of pulse periods found in accreting pulsar 
systems (at least in the SMC). There is obviously a strong 
preference to forming systems with spin periods of <100s. 
Since the Corbet Diagram (Corbet et al., 2008) provides 
strong evidence that the equilibrium spin period is driven 
by the orbital period, then this spin period result probably 
reflects the relative difficulty in long orbital period systems 
surviving the supernova explosion that produces the neutron 
star. This is supported by the binary period distribution seen 
in the lower panel of Figure |3] (but with all the observational 
caveats mentioned above). It is also worth noting that there 
is no evidence from the Corbet diagram that SMC sources 
show any difference from their counterparts in the Milky 
Way (Corbet et al., 2008). 



in the SMC. This homogenous group of objects provides an 
excellent test bed for accretion theory and, in general, the 
results are very supportive of current models. Since current 
models have all been developed based upon Milky Way sys- 
tems, this work strongly supports the conclusion that there 
is little, or no difference between HMXB behaviour in the 
SMC and the Galaxy. 
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4 CONCLUSIONS 

This work presents results from -^10 years of X-ray mon- 
itoring of the pulse period histories of 24 HMXB systems 
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